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^ We have just entered an era of precision measurements for S^-decay observ- 

^ ables. A characteristic feature of the S^-meson system is B^-B^ mixing, 

which exhibits a sizable decay width difference. The latter feature leads to a 
O subtle complication for the extraction of branching ratios of Bg decays from 

^ untagged data samples, leading to systematic biases as large as O{10%) that 

depend on the dynamics of the considered decay. We point out that this 
•'-j effect can only be corrected for using information from a time-dependent 

rS analysis and suggest the use of the effective B^ decay lifetime, which can 

^ already be extracted from the untagged data sample, for this purpose. We 

also address several experimental issues that can play a role in the extrac- 
tion of effective lifetimes at a hadron collider, and advocate the use of the 
Bg branching ratios, as presented in this note, for consistent comparisons of 
theoretical calculations and experimental measurements in particle listings. 
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We have just entered an era of precision measurements for _Bs-decay observables. A characteristic 
feature of the Bs-meson system is Bg-B'^ mixing, which exhibits a sizable decay width difference. 
The latter feature leads to a subtle complication for the extraction of branching ratios of Bs decays 
from untagged data samples, leading to systematic biases as large as 0(10%) that depend on the 
dynamics of the considered decay. We point out that this effect can only be corrected for using 
information from a time-dependent analysis and suggest the use of the effective Bs decay lifetime, 
which can already be extracted from the untagged data sample, for this purpose. We also address 
several experimental issues that can play a role in the extraction of effective lifetimes at a hadron 
collider, and advocate the use of the Bs branching ratios, as presented in this note, for consistent 
comparisons of theoretical calculations and experimental measurements in particle listings. 
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I. INTRODUCTION 

Weak decays of Bg mesons encode valuable informa- 
tion for the exploration of the Standard Model (SM). The 
simplest observables are branching ratios, which give the 
probability of the considered decay to occur. Measure- 
ments of Bs branching ratios at hadron colliders, such as 
Fermilab's Tevatron and CERN's Large Hadron Collider 
(LHC), would require knowledge of the Bg production 
cross-section, which presently makes absolute branch- 
ing ratio measurements impossible. Hence experimental 
control channels and the ratio of the fs/ fu,d fragmenta- 
tion functions, describing the probability that a b quark 
hadronizes as a Bq meson [T] , are required for the conver- 
sion of the observed number of decays into the branch- 
ing ratio. At e+e^ B factories operated at the T{5S) 
resonance, the total number of produced Bg mesons is 
measured separately and subsequently also allows for the 
extraction of the Bg branching ratio from the data [5] . 

A key feature of the Bg mesons is B^-B^ mixing, 
which leads to quantum-mechanical, time-dependent os- 
cillations between the B^ and B^ states. In contrast to 
the Bd system, the Bg mesons exhibit a sizable differ- 
ence between the decay widths of the light and heavy 
mass eigenstates, r^"'' and F^'', respectively [3]. Cur- 
rently the most precise measurement is extracted from 
the — J/'ipi' channel by the LHCb collaboration [i]: 
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0.088 ±0.014; 



(1) 



T^^>]/2 = (0.6580 ± 0.0085) ps-i is 
the inverse of the Bg mean lifetime tb^ ■ 

In view of the sizable decay width difference, Eq. ([T]), 
special care has to be taken when dealing with the con- 
cept of a branching ratio. We shall clarify this issue and 
give an expression, allowing us to convert the experimen- 
tally measured Bg branching ratio into the corresponding 



"theoretical" branching ratio. The latter is not affected 
by B^-B^ mixing and encodes the information for the 
comparison with branching ratios of B^ decays, where 
the relative decay width difference at the 10"'^ level [3] 
can be neglected, or branching ratios of i3+ modes. 

The difference between these two branching ratio con- 
cepts involves yg and is specific for the considered Bg 
decay, thereby involving non-perturbative parameters. 
However, measuring the effective lifetime of the consid- 
ered Bg decay, the effect can be included in a clean way. 

In experimental analyses, this subtle effect has so far 
been neglected or only been partially addressed; exam- 
ples are the branching ratio measurements of the Bg — >■ 
K+K- Bg ^ J/i^foi980) Bg ^ J/^Kg H, 
Bg D+D- and DjTr+ ^ decays by the 

LHCb, CDF, D0 and Belle collaborations. 



II. EXPERIMENT VERSUS THEORY 

What complicates the concept of a Bg branching ratio 
is the fact that the untagged decay rate is the sum of two 
exponentials [10] : 

{TiBgit) ^ /)) ^ TiB'^git) -^f) + rmt) -> /) 
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corresponding to two mass eigenstates with different life- 
times. Using Eq. Q, we write 
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is a final-state dependent observable. 

In experiment it is common practice to extract a 
branching ratio from the total event yield, ignoring in- 
formation on the particles' lifetime. The "experimental" 
branching ratio can thus be defined as follows [TU] : 
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r{Bs{t) ^ f)) dt (5) 



rL + r 



Note that this quantity is the average of the branching 
ratios for the heavy and light mass eigenstates. 

On the other hand, what is generally calculated the- 
oretically are CP-averaged decay rates in the flavor- 
eigenstate basis, i.e. 

(r(i?.(t) ^ - r(i?° ^ /) + r(B," ^ /). (6) 

This leads to the following definition of the "theoretical" 
branching ratio: 
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By considering i = 0, the effect of B^^-Bl mixing is 
"switched off'. The advantage of this Bg branching 
ratio definition, which has been used, for instance in 
Refs. [HHH], is that it allows a straightforward compar- 
ison with branching ratios of B^ or B^ mesons by means 
of the SU{2i) flavor symmetry of strong interactions. 

The experimentally measurable branching ratio, 
Eq. ^, can be converted into the "theoretical" branch- 
ing ratio defined by Eq. Q through 



BR(i3,^/)^ 



(8) 

In the case of ys — 0, the theoretical and experimental 
branching ratio definitions are equal. 

Inspection of Eq. ([8| reveals that ys and -4^p are re- 
quired for the translation of the experimental branching 
ratios into their theoretical counterparts. Ideally, the lat- 
ter quantities should eventually be used in particle com- 
pilations, in our opinion. 

The decay width parameter y^ is universal and has 
already been measured, as summarized in Eq. ([I]). In 
Fig. [l] we illustrate Eq. Q for a variety of values of 
and observe that differences between BR [Bs — > /)thco 
and BR [B^ /)exp large as 0(10%) may arise. 

The simplest situation corresponds to flavor-specific 
(FS) decays such as i?° — > D~it^ , where = and 
the correction factor is simply given by 1 — y^. 

However, if both the J5° and the B? mesons can de- 



cay into the final state /, the observable -4^p is more 
involved and depends, in general, on non-perturbative 
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FIG. 1: Illustration of Eq. for various values o/^^p. We 
also show the current LHCb measurement of ys 



hadronic parameters, CP-violating weak decay phases, 
and the B^-B^ mixing phase cjis ■ Assuming the SM struc- 
ture for the decay amplitudes and using the SU{Z) fiavor 
symmetry to determine the hadronic parameters from 
relations to Bd decays, theoretical analyses of ^^p were 
performed for the final states J/tpcj) [E], K^K~ [13], 
J/V^/o(980) mi, J/ijKs Hg and DfD- IF]. 



III. USING LIFETIME INFORMATION 



The simplest possibility for implementing Eq. Q is to 
use theoretical information about the -4£p observables. 
However, this input can be avoided once time information 
of the untagged Bg decay data sample becomes available. 
Then the effective lifetime of the Bg ^ f decay can be 
determined, which is theoretically defined as the time 
expectation value of the untagged rate |17) : 



^ S^t{T{Bg{t)^f))dt 
j^{T{Bg{t)^ f)) dt 



l + 2Ai^ys+yl 



(9) 



The advantage of Tf is that it allows an efficient extrac- 
tion of the product of y^£p and j/j. Using the effective 
lifetime, Eq. Q can be expressed as 



BR{Bg^f) 
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(10) 

Note that on the right-hand side of this equation only 
measurable quantities appear and that the decay width 
difference yg enters at second order. The measurement of 
effective lifetimes is hence not only an interesting topic 
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BR{Bs /)cxp 
(measured) 



^ir(SM) 



(1.29t°;«) X 10-* m 
(3.5 ±0.8) X 10"^ [7] 
(3.01 ±0.34) X 10"^ [9] 
(3.5 ±0.7) X 10"^ [IH] 

(1.041^2^) X 10-^ m 



From Eg. ((8| 



From Eq. ( 10 I 



J /i>f 0(980) 

K+K- 
D+D- 



0.9984 ±0.0021 [H] 
0.84 ±0.17 dH 
(exact) 
-0.972 ±0.012 [13] 
-0.995 ±0.013 [IS] 



0.912 ±0.014 
0.924 ±0.018 
0.992 ± 0.003 
1.085 ±0.014 
1.088 ±0.014 



0.890 ±0.082 j6j 

N/A 

N/A 
1.042 ±0.033 dn] 

N/A 



TABLE I; Factors for converting BR{Bs /)^^p (see ^) into BR (B^ ./■)thco (^ee Q) by 

means of Eq. ([sjl with theoretical 



estimates for .4£p. Whenever effective hfetime information is available, the corrections are also calculated using Eq. (10 1 



for obtaining constraints on the B^-B'^ mixing parame- 
ters [T^ , but an integral part of the determination of the 
"theoretical" B^ branching ratios from the data. 

In Table |lj we hst the correction factors for converting 
the experimentaUy measured branching ratios into the 
theoretical branching ratios for various decays. Here we 
have used theoretical information for .A^p and Eq. (l8), or 
- if available - the effective decay lifetimes and Eq. ( 10| . 

The rare decay B'^ — >■ fj,^fj,~, which is very sensitive to 
New Physics [2D], is also affected by AFg. In Ref. [21,, we 
give a detailed discussion of this key Bg decay, showing 
that the heUcities of the muons need not be measured to 
deal with this problem, and that AF^ actually offers a 
new window for New Physics in i?^ 



IV. 



VV DECAYS 



Another application is given by Bs transitions into two 
vector mesons, such as Bs -J> J/^(j> |22J, -> x*'^K*'^ 
[55] and B^ D*'^D*~ [E]. Here an angular analy- 
sis of the decay products of the vector mesons has to 
be performed to disentangle the CP-even and CP-odd 
final states, which affects the branching fraction deter- 
mination in a subtle way, as recognized in Refs. [23J . 
Using linear polarization states 0, || with CP eigenvalue 
rjk = +1 and ± with CP eigenvalue rjk = —1 [IS], the 
generalization of Eq. (l8|) is given by 
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As discussed in Ref. 
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ivv,k — ^- uiDuuDocu 111 _Lici.. [T7] , assummg 
the SM structure for the decay amplitudes, we can write 
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(12) 



where Cyv.k describes direct CP violation, (jjg is the 
B^-B^ mixing phase, and A(j)vv,k is a non-perturbative 
hadronic phase shift. The expressions given in Ref. [23] 
for the Bs — K*'^K*'^ decay take the leading order ef- 
fect of Us into account, and assume (f)s — and negligible 
hadronic corrections. 



The generalization of Eq. (|10|) is given by 

^^thco — ^^cxp / , 



k=0.\\,± 



2 - (1 - y's) 



2\ ^k 



rcxp 

JVV.k' 



(13) 

and does not require knowledge of the A^y' observables. 



V. EXPERIMENTAL ASPECTS 

Additional subtleties arise in the experimental deter- 
mination of Bs branching ratios and effective lifetimes, in 
particular at a hadron collider environment where many 
final-state particles are produced in the fragmentation. 

Separating Bs signal decays from the background typi- 
cally involves selection criteria that use the flight distance 
of the Bs meson or the impact parameter of its decay 
products, leading to a decay-time dependent efficiency. 
By rejecting short-living Bs meson candidates, the rela- 
tive amounts of Bs^h and i?s,H mesons in the remaining 
data sample are altered, resulting in a biased result for 
the branching ratio determination. The extrapolation of 
the event yield to full acceptance is usually obtained from 
simulation, but this requires a priori assumptions of the 
values for ys and .A^p. For example, the dependence of 

the branching fraction correction on the value ^^p can 
be several percent if only decay times greater than 0.5 ps 
are considered. This systematic uncertainty is avoided 
by tuning the simulation using the measured value of the 
effective lifetime. 

Furthermore, the presence of remaining background 
events with a different observed decay time distribution 
as the signal, implies that it is experimentally unprac- 
tical to determine the time expectation value Tf of the 
untagged rate as given in Eq. ([9|. Instead, the effective 
lifetime is commonly extracted by fitting a single expo- 
nential function to the untagged rate [S] HH US] , which 
in general is described by two exponentials (see Eq. ([2|). 
In Appendix|A]we demonstrate that such a fitting proce- 
dure leads to an unbiased determination of the effective 
lifetime in the case of a log likelihood fit and to a small 
bias for a minimization procedure. 
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VI. CONCLUSIONS 

The established width difference of the Bg mesons com- 
plicates the extraction of branching ratio information 
from the experimental data, leading to biases at the 10% 
level that depend on the specific final state. On the one 
hand, these effects can be included through theoretical 
considerations and phenomenological analyses. On the 
other hand, it is also possible to take them into account 
through the measurement of the effective i?s — / decay 
lifetimes, which is the preferred avenue. So far, these ef- 
fects have not, or only partially, been included and we 
advocate to use the converted branching ratios for com- 
parisons with theoretical calculations in particle listings. 



Appendix A: Effective Lifetime Fits 

An effective lifetime for a Bg decay channel is obtained 
in practice by fitting a single exponential function to its 
untagged rate. As an untagged rate is in general de- 
scribed by two exponentials, corresponding to two mass- 
eigenstates with different lifetimes, the single exponential 
fit is an approximation. 

In order to find analytic expressions for the fitted ef- 
fective lifetime Toff, we let the untagged rate be the true 
Probability Distribution Function (PDF), and the single 
exponent function the fitted PDF, such that 



^ A{t) (Pft)) 

/trucW- ^oo^^^^ {T{t))dt' 



/fit(i; Tcff) 



A{t) e 



J^-'A{t)e'^/r^«dt' 



(Al) 
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where A{t) is an acceptance efficiency function. The like- 
lihood or function for the fit in question is then built 
using the above PDFs, and maximised or minimised, 
respectively, in the limit of infinitcsimally spaced bins. 



Specifically, for n events we minimise the functions: 

/>oo 

- log L{tcs) = - n / dt /truc(t) log [/fit(i; Tcff )] , 
Jo 

(A3) 

2f X f°°,,[fti-uc{t)-f{it{t;Tcs)f , , 

X (TcB) = n / dt — r , (A4) 

for a maximum likelihood and a least squares fit, respec- 
tively. In a modified least squares fit, where data is used 
to estimate the error, the denominator in the inte- 
grand should be replaced by /truo(i)- For the maximum 
likelihood fit, taking the infinitesimal bin limit is equiv- 
alent to an unbinned fit. 

The effective lifetime Tcs resulting from these fits is 
then given implicitly by the formula: 



/„°°te-*/"°"A(t) dt ^ J^tg{t;T,g)A{t) dt 
j;^e-t/-^«A{t)dt ~ j;;°g{t;T,s)A{t)dt' 



(A5) 



where 



l{t] Toff) 



im) 

(r(i))2e*/^''« 
(r(i))-ie-2t/^cff 



maximum likelihood 
least squares 
modified least squares. 



The effective lifetime definition given in ([9| is repro- 
duced for the untagged rate given in (H if we assume a 
trivial acceptance function, A{t) — 1, and apply a max- 
imum likelihood fit. For non-zero values of yg, the least 
squares fits give different analytic expressions for the ef- 
fective lifetime. Fortunately, for the current experimental 
range of ys, the differences are of the order 0.1%. 
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